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a solution to the Hermitian Yang-Mills equations. We then extend our results by 
investigating the behavior of non slope-stable bundles. In a variety of examples, 
it is shown that the failure of these bundles to satisfy the Hermitian Yang-Mills 
equations, including field-strength singularities, can be accurately reproduced 
numerically. These results make it possible to numerically determine whether 
or not a vector bundle is slope-stable, thus providing an important new tool in 
the exploration of heterotic vacua. 



Email: andlara@physics.upenn.edu, vbraun@stp.dias.ie, rlk@vt.edu, 

ovrutOelcapitan . hep . upenn . edu 



Contents 



1 Introduction 

1.1 The Generalized Donaldson Algorithm and Heterotic Compactifications 

1.2 The Hermitian Yang-Mills Equations and Hermitian- Einstein Bundle 
Metrics 

1.3 Slope Stability 

2 Computing the Calabi-Yau Metric 

2.1 Donaldson's Algorithm 

2.2 Ricci Curvature and Error Measures 

2.3 The Quartic K3 

2.4 The Quintic Threefold 

3 Hermitian Yang-Mills Connections 

3.1 Generalizing Donaldson's Algorithm 

3.2 Untwisting 

3.3 Measuring the Error 

3.4 An Example 

3.5 Subtracting the Trace 

3.6 Integrated Error Measure 

4 Stable vs. Unstable 

4.1 Taxonomy of Slope-Stability 

4.2 The Harder-Narasimhan Filtration 

4.3 Examples of Semi-Stable and Unstable Bundles on the Quartic K3 . . 

4.3.1 A Semi-Stable Bundle 

4.3.2 An Unstable Bundle with a Filtration by Line Bundles 

4.3.3 An Unstable Bundle with Sheaf Filtration 

4.4 Eigenvalues Along a Geodesic 

4.4.1 Iterating the T-Operator 

4.4.2 Degree of Twisting 

4.5 Rank Three Bundles on the Quintic 

4.5.1 A stable bundle 

4.5.2 A semi-stable bundle 

4.5.3 Unstable bundles 

5 Conclusions and future work 
Bibliography 



1 



1 Introduction 



1.1 The Generalized Donaldson Algorithm and Heterotic Com- 
pact ificat ions 

A central goal of string theory is to produce low-energy theories with the symmetries, 
spectrum and properties of elementary particle physics. In the search for realistic 
vacua, compactification of the heterotic string [H [2] and M-theory [31 HJ E] on smooth 
Calabi-Yau threefolds with holomorphic vector bundles [61 [7J El El [TUl [HJ EE2] has 
played an important role [HI HH [15l [16j [17] . 

Heterotic compactifications possess a number of phenomenologically desirable fea- 
tures, including realistic gauge groups and particle spectra [HI [191 EQJ EH [22} [231 
EH [25] as well as gauge coupling unification [261 [27] . However, certain fundamen- 
tal quantities, such as the Yukawa couplings, are difficult to compute directly [28l 
|29"1 [30] [3T] [32] 133] . For example, the "physical" Yukawa couplings depend on both 
the Yukawa coefficients in the superpotential and the explicit form of the Kahler po- 
tential. Both quantities are determined by the detailed structure of the underlying 
geometry-specifically, the metric on the Calabi-Yau threefold and the connection on 
the slope-stable holomorphic vector bundle-about which, generically, little is known. 
It follows that the Yukawa couplings in the four- dimensional effective theory have not 
been explicitly computed, except in very special cases where sophisticated tools from 
algebraic geometry or topological string theory are available. As a result, one can 
rarely do better than the qualitative statement that such coupling coefficents either 
vanish or are "of order one" . 

To fully specify the geometry and, hence, to be able to compute the couplings 
in the effective theory, one must determine the Ricci-flat metric g on the Calabi- 
Yau threefold X and the supersymmetric gauge connection A on the slope-stable 
holomorphic vector bundle V [31] . Yau's theorem [32] gives an existence proof for the 
Ricci-flat metric and the work of Donaldson-Uhlenbeck-Yau [361 EZ] provides us with 
a class of holomorphic vector bundles, called slope-stable, that are consistent with 
supersymmetry. However, in general the form of g and A are not known analytically. 
The technical challenge of determining these quantities has made systematic studies 
of realistic heterotic vacua difficult to achieve. 

In recent years, a number of new numerical approaches to these old problems have 
been presented [381 EH SQl SH g3l SH US SSI gTJ. Thanks to the development of 
powerful new algorithms and increased computer speed, it is now possible to compute 
Ricci-flat metrics and to directly solve the Hermitian Yang-Mills equations for the 
gauge connection. With this data in hand, all quantities in the four-dimensional 
effective theory, such as the correctly normalized zero modes, all coefficients in the 
superpotentials and the Kahler potential, can, in principle, be computed. In this 
paper, we take several important new steps towards this goal. The paper is structured 
as follows. 
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In Section El we outline the general form of Donaldson's algorithm for computing 
the Ricci-flat metric [381 EHl HO] on a Calabi-Yau manifold. Specifically, using the 
numerical implementations developed in [5HI SS] and [SHI S3], we calculate Ricci-flat 
metrics on the K3 surface and the Quintic threefold. To evaluate the accuracy of our 
approximation scheme, that is, how close our approximate metric is to the Ricci-flat 
solution, several measures of the metric-error are introduced. In Section [HI we turn 
our attention to recent generalizations of Donaldson's algorithm due to Yau, Wang 
and Douglas et. al. [121 E2]- These make it possible to compute Hermitian-Einstein 
metrics on holomorphic vector bundles over a Calabi-Yau manifold and, hence, to 
solve for the gauge connection satisfying the conditions for M = 1 supersymmetry. 

Given a rank n holomorphic vector bundle V on a Calabi-Yau manifold X, the 
central idea of the algorithm rests on establishing an embedding of X into a Grass- 
manian G(n,Nk — 1) via the Nk sections of V ® C® k , where £ is some ample line 
bundle. Using this embedding, a family of simple Kahler metrics on G(n, Nk — 1) can 
be pulled-back to the Calabi-Yau manifold and used to define a sequence of Hermi- 
tian fiber metrics converging to a Hermitian-Einstein metric on V. In [^2l [50] . Wang 
proposed a generalization of Donaldson's "T-operator" for the embedding described 
above. It was shown that the iteration of this T-operator converges to a fixed point for 
each value of k if and only if the bundle V is Gieseker- stable. Furthermore, in the case 
that V is slope-stable, there exists a scheme of approximate fiber metrics on V which 
converges to the Hermitian-Einstein metric in the limit that k — > oo. Wang's results 
are reviewed in Section 13.11 We then present a new algorithm to explicitly construct 
the Hermitian-Einstein fiber metric and, hence, the supersymmetric gauge connection 
A on V. An interesting technical aspect arises in the process of extracting the con- 
nection on V from the connection produced by the algorithm on the "twisted" bundle 
V®£® k . In Section [3721 by studying all possible methods of "untwisting" the bundle, 
we establish the computationally optimal choice and present a general procedure for 
the computation of the connection on V itself. Following Wang's theorem (Theorem 
H]in Section I3~f]) . we explicitly compute the connection on several slope-stable vector 
bundles, beginning in Section \3.3[ and develop an accurate error measure, Tk H , to 
quantify the accuracy of numerical approximations. 

In Section HJ we move beyond Wang's theorem and explore the behavior of non- 
stable bundles under the generalized Donaldson algorithm. Since Wang's theorem does 
not apply in this case, there is, a priori, no reason to expect any specific behavior 
or systematic results from the algorithm. We show, however, that the output of 
Donaldson's algorithm, when applied to properly semi-stable or unstable bundles, can 
in be understood in terms of their Harder- Narasimhan filtrations. Furthermore, it is 
shown in detail that the failure of these bundles to satisfy the Hermitian Yang-Mills 
equations, ranging from non-trivial constants to field-strength singularities, can be 
accurately reproduced numerically. These observations demonstrate that Donaldson's 
algorithm provides an important new tool in the study of the slope-stability properties 
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of vector bundles. In an arbitrary heterotic compactification, one of the most difficult 
obstacles to overcome is deciding whether or not a given vector bundle is slope-stable. 
This determination has hithertofore required analytically determining all subsheaves 
of the bundle and calculating their slopes. This, at best, is a difficult thing to do, 
and often cannot be carried out. In this paper, we show that the stability properties 
of a vector bundle, that is, whether it is stable, semi-stable or unstable, can always be 
determined by numerical calculation using the generalized Donaldson's algorithm. This 
result provides a new tool in the exploration of supersymmetric heterotic vacua. We 
illustrate these procedures by looking at a variety of stable, semi-stable, and unstable 
bundles on the K3 twofold and the Quintic threefold. The bundles we choose as 
examples are defined via the monad construction [5H [52} [53], [25], a class known to 
contain bundles corresponding to physically realistic heterotic compactifications. 

In the following two subections, we set the stage for the entire discussion by defin- 
ing the Hermitian Yang-Mills equations and their relationship to the slope-stabilty of 
vector bundles [54] . 

1.2 The Hermitian Yang-Mills Equations and Hermit ian-Ein- 
stein Bundle Metrics 

An M = 1 supersymmetric heterotic string compactification is specified by a com- 
plex three-dimensional Calabi-Yau manifold, X, and a holomorphic vector bundle, V, 
with structure group K C E$ defined over X. The gauge connection, A, on V with 
associated field strength, F, must satisfy the well-known equations [M] 

F ij = F l3 = 0, g fj F fj = (1.1) 

where i, j = 1,2,3 run over the holomorphic indices of the Calabi-Yau threefold. These 
equations arise as the vanishing condition on the variation of the 10- dimensional gaug- 
ino and, hence, are required to preserve M — 1 supersymmetry. The first two condi- 
tions, Fij = Fjj = 0, are simply the constraint that the vector bundle be holomorphic. 
On a holomorphic vector bundle, that is, one with holomorphic transition functions, 
one can always choose a connection with a purely (1,1) field strength. Hence, the van- 
ishing of the (2, 0) and (0, 2) components is satisfied. The last condition, g^Fq = 0, 
however, is not so easily solved. Equations (11.11) are a special case of the Hermitian 
Yang-Mills equations. More generally, one has 

g q F fj = c-\ (1.2) 

where c is a constant determined by the first Chern class of the bundle V. For a generic 
holomorphic, rank n vector bundle, the structure group is realized in the fundamental 
representation of U(n). Then Ci(V) ^ and, therefore, c ^ 0. However, for strucure 
group K C E s , it follows that ci(V) must vanish; hence, the zero trace in (11. ip . The 
numerical algorithms in this paper will require us to find solutions to the general 
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equation ( 11.21) . The results will then be specializing to the case c = 0, that is, ( II. ip . 
relevant for heterotic compactifications. 

A solution to ( II. 2p is equivalent to the bundle V carrying a particular Hermitian 
structure. An Hermitian structure (or Hermitian fiber metric), G, on V is an Her- 
mitian scalar product G x on each fiber V(x) which depends differentiably on x. The 
pair (V, G) is often referred to as an Hermitian vector bundle. For a given choice of 
frame, e a (x), for the bundle, one can define the covariant derivative as 

D(v a e a ) = (dv a )e a + v a A b a e b , (1.3) 

where a,b = 1, . . .n take values in the structure group K C U(n). The Hermitian 
structure defines an inner product as 

(e a ,e b ) = G- ab , G = G\ (1.4) 

The condition that the connection be compatible with the metric, 

d(u, v ) = (Du, v ) + (u, Dv) , (1.5) 

gives 

dG- ab = G- ac A c b (1.6) 

and, hence, 

A = G- l dG . (1.7) 

Since one can always make a gauge choice to set A = for a holomorphic bundle, we 
can re-phrase the Hermitian Yang-Mills equation ( 11.21) on i^ 1,1 ) as a condition on the 
metric, 

c • 1 = <7% = = g^ i dj(G~ 1 d i G) . (1.8) 

A metric G on the fiber of satisfying the above equation is called an Hermitian- 
Einstein metric. 

Finally, note that given a fiber metric G, one can define an inner product on the 
space of global sections, where a — 1, . . . h°(X, V), of V as 

< s a \s p >= [ {s b pG b - a sl) dVol . (1.9) 
Jx 

This will be of use to us in subsequent sections. Before attempting to explicitly 
construct a connection A on V satisfying (\1.2\\ (or, equivalently, an Hermitian-Einstein 
bundle metric satisfying (I1.8P ). we will first rephrase these conditions in terms of 
algebraic geometry and recall a central theorem. 
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1.3 Slope Stability 



Having presented the Hermitian Yang-Mills equations and defined Hermitian-Einstein 
bundle metrics, we now introduce one more important mathematical notion. As stated 
above, the Hermitian Yang-Mills equation g^F^ = c- 1 is notoriously difficult to solve 
for the case of non-Abelian structure groups. However, on Calabi-Yau manifolds there 
exists a powerful method for transforming this equation into a problem in algebraic ge- 
ometry. For Kahler manifolds, the Donaldson-Uhlenbeck-Yau theorem [36j [37] states 
that on each poly-stable holomorphic vector bundle V there exists a unique connec- 
tion satisfying the Hermitian Yang-Mills equation (jl.2p . Thus, to verify that a vector 
bundle is consistent with supersymmetry, one need only verify that it possesses the 
property of poly-stability. 

The notion of slope-stability (also referred to as Mumford-Takemoto stability) of 
a bundle (or coherent sheaf), J 7 , over a Kahler manifold, X, is defined by means of a 
real number called the slope: 

H{T) = I Cl (F)Au d -\ (1.10) 

rk(.F) Jx 

where d is the dimension of the Kahler manifold. Here, u is the Kahler form on X, 
while rk(J r ) and ci(jF) are the rank and the first Chern class of J 7 respectively. A bun- 
dle V is called stable (semi-stable) if for all sub-sheaves J 7 C V with < rk(J r ) < rk(V) 
the slope satisfies 

^)<fi(V) (//(J") < n(V)) . (1.11) 

A bundle is poly-stable if it can be decomposed into a direct sum of stable bundles 
which all have the same slope. That is, 

v = 0v n , pcvy = a*(V) . (i.i2) 

n 

Thus, stability is a special case of poly-stability which is in turn a subset of semi- 
stability. As a series of implications: stable =>- poly-stable =>■ semi-stable. It is 
important to note that the converse to these statements do not hold. For example, 
not every semi-stable bundle is poly-stable. Finally, observe that the slope is exactly 
the constant, c, that appeared in (II .2|) . Written in terms of the slope of V, the general 
Hermitian Yang-Mills condition is 

g* 3 Fy=ii(V)-l nxn , (1.13) 

where n is the rank of V and ^(V) is defined in (11 .101) . 

It should be noted that the condition of slope-stability on a Calabi-Yau manifold 
depends on all moduli of the heterotic compactification. To be specific, consider a 
Calabi-Yau threefold. Here, the moduli are the Kahler moduli (/i 1 ' 1 (X)), the com- 
plex structure moduli (h 2,1 (X)) and the vector bundle moduli (h 1 (End(V))) . The 
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dependence on Kahler moduli is explicit in ( 11. ip and referred to as a choice of "po- 
larization". One can expand the Kahler form u in (II. 1 f) as u = t r u r , where u r 
form a basis of (1, l)-forms and t r are the real parts of the Kahler moduli. Written 
in terms of the triple intersection numbers, d rst , of the threefold, the slope is simply 
A*(T) = '^yjd rs tCi(V) r t s t t where r, s = 1, . . . The other moduli enter through 

the notion of a subsheaf J G V. In general, whether or not there exists an injective 
sheaf morphism : T — > V depends on both the complex structure moduli and the 
vector bundle moduli. 

Thus, finding a solution to the Hermitian Yang-Mills equations (that is, deter- 
mining whether the bundle is slope-stable) is a question that must be asked after 
selecting a particular point in moduli space. As we will see, some regions of moduli 
space may admit a solution while others will not. This moduli dependence can lead 
to a variety of interesting physical consequences [551 EH1 EU [58], which we will explore 
in later sections. For most of this work, however, we will be interested in studying 
the stability properties of bundles at "generic" points in their moduli space. 

Before investigating solutions to the Hermitian Yang-Mills equations, we must 
first introduce an algorithmic approach to determining Hermitian metrics (both of 
manifolds and bundles). In the following section, we discuss Donaldson's algorithm 
for numerically approximating the Ricci-flat metric on a Calabi-Yau manifold. 

2 Computing the Calabi-Yau Metric 
2.1 Donaldson's Algorithm 

Many of the challenges associated with string compactifications on a Calabi-Yau three- 
fold, X, arise from the difficulty in determining the explicit geometry. The simplest 
M = 1 supersymmetric vacuum solutions require a Ricci-flat metric, g^, on X. While 
Yau's theorem [55] ensures that such a metric exists, no analytic solutions have been 
found. Recently, however, it has become possible to find numerical solutions using an 
algorithm developed by Donaldson [38j |39j 00] . This algorithm has been implemented 
numerically and extended in [59 | B8 | B9 | |50 | H3 | H6] . For other numerical computations 
of Kahler metrics, see [lU EDI US]- In this paper, we propose to investigate solutions 
to the Hermitian Yang-Mills equation g^Fq = 0. Since the Ricci-flat metric, g^, is 
required for such a solution, in this section we provide a brief review of Donaldson's 
algorithm for approximating such metrics on Calabi-Yau manifolds. 

The starting point of Donaldson's algorithm is the observation that there exists 
a natural metric on P n . Denoting the n + 1 homogeneous coordinates by Zi, this 
Fubini-Study metric is given by gpsfj = ^didjK FS , where 




(2.1) 



7 



and h 1 -* is any hermitian, non-singular matrix. The Fubini-Study metric is usually 
defined with = but we will refer to ( 12. ip as a Fubini-Study Kahler potential 
for any hermitian matrix h % K The Fubini-Study metric can be used to induce a 
metric on any subvariety of P n . In particular, it will induce a metric on any Calabi- 
Yau threefold X, since it is always possible to embed IcP" for some n. However, 
such a metric will not, in general, be Ricci-flat. To obtain the Ricci-flat metric, 
Donaldson's algorithm requires one to use a generalized version of ( 12. ip with enough 
free parameters to provide a more versatile induced metric on X. The algorithm then 
guarantees that this generalized metric will converge, in some specified limit, to the 
Ricci-flat metric. 

An obvious generalization of ( 12. ip is to replace the degree one polynomials with 
polynomials of higher degree. For example, 

K=^ E h tl - lk]l - 3k z H ...z tk z ]l ...z ]k (2.2) 

h—ikh—Ok 

where h %1 "' %k ^ 1 "^ k is hermitian. This new Kahler potential now has (n + l) 2k real pa- 
rameters. One can write the above generalization in a more systematic way by noting 
that such Kahler potentials can naturally be obtained by considering holomorphic 
line bundles on the Calabi-Yau manifold X itself and using the Kodaira embedding 
theorem. Let L be a holomorphic line bundle over X with n\ = h°(X,L) global 
sections. Consider a twisting of the line bundle, C k = JZ® k . Then, choosing a basis 
for the space of sections, s a G H°(X,L k ) where < a < — 1, allows one to define 
a map from X to P nfc_1 given by 

i k : X -> F*" 1 , (x , • • • , x 2 ) h-> [s Q (x) : . . . : s nfc _i(z)] , (2.3) 

where Xi are holomorphic coordinates on the Calabi-Yau threefold. The Kodaira 
embedding theorem [61] states that if C is ample, (12.31) will define an embedding of 
X C P"*" 1 for all C k with k > k for some k . 

In terms of this natural embedding line bundle, we can view the generalized Kahler 
potential (12. 2p restricted to X simply as 

n/j — 1 

Kh > k = E ^ a « 5 * = ^Mlk ■ (2-4) 

o,/3=0 

Geometrically, ( 12. 4p defines an hermitian metric on the line bundle L k itself. It 
provides a natural inner product on the space of global sections 

= (s p \s a ) = v] Hk (Y , [ ^| dVoW , (2.5) 

VOlcY{X) J x \ \S\\ h 

where 

dVolcY = A n (2.6) 
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and Q is the holomorphic (3,0) volume form on X . Note that ( 12.5ft depends non- 
linearly on h. It was shown by Tian jH] that metrics on ample line bundles, such 
as (12.41) . can provide a useful "basis" of Kahler metrics on X. That is, any algebraic 
Kahler potential can either be written in the form ( 12.41) or obtained as the limit point 
of a sequence of such potentials. Specifically, one has the following theorem. 

Theorem 1 (Tian). Let {s a } be a basis for H°(X,L k ) for some ample line bundle 
L. Then the space of all "algebraic" Kahler potentials, 

1 ^fc -1 

KIT _ 

a,/3=0 

where k e 7L, is dense in the space of Kahler potentials. 

With this observation, we return to the goal of finding the Ricci-flat metric on 
X. Consider X as a projective manifold, with an ample line bundle £. Note that, 
in general, the matrices h al3 and M a p in (12. 5p are completely unrelated. However, 
for special metrics, they may coincide. The metric h on the line bundle L is called 
"balanced" if 

(M^)" 1 = h a K (2.8) 

In this case, one can find an "orthonormal" basis of sections for which M a p = 5 a p and 
h a/3 = 5 al3 . Many theorems have been proven about balanced metrics. Of particular 
interest here is their curvature properties [38| |39| HQ]. 

Theorem 2 (Donaldson). For each k > 1, the balanced metric, h, on L k exists and 
is unique. As k — >■ oo, the sequence of metrics 

a,j3=0 

on X converges to the unique Ricci-flat metric for the given Kahler class and complex 
structure. 

The core of Donaldson's algorithm is then simply the task of finding the balanced 
metric for each k. To this end, Donaldson defines the "T-operator" 

T{hU = ^TJxS I v S S - dVoW (2 - 10) 

Vol C y(A) J x ^gh^SjSs 

which, for a metric h, computes the matrix T(h). Given a fixed point of this operator, 
that is, a metric h for which T(h) = h, one can always perform a change of basis 
s — > h~ l l 2 s to bring h to the unit matrix, producing a balanced embedding. To find 
this fixed point, simply iterate ( I2.10p as follows. 



9 



Theorem 3 (Donaldson). For any initial metric ho (and basis s a of global sections 
of L k ), the sequence 

h n+1 = (TiK))" 1 (2.11) 
converges to the balanced metric as n — > oo. 

Happily, in practice, very few (< 10) iterations are needed to approach the fixed point. 
Henceforth, we will also refer to g^ 1 in (12.91) . the approximating metric for fixed k, as 
a balanced metric. 

To find the balanced metric at each step k, one must be able to integrate over the 
Calabi-Yau threefold. An explicit numerical integration scheme for this purpose was 
given in [421 EH] and is used to compute the Calabi-Yau metrics in this paper. We 
refer the reader to these references for details of this and other technical aspects of 
the computer implementation. To summarize Donaldson's algorithm: 

1. Choose an ample line bundle L and a degree k (that is, a twisting of the line 
bundle C k ) at which to compute the balanced metric which will approximate 
the Calabi-Yau metric. 

2. Calculate a basis {sq}™^ 1 for H°(X,C k ) at the chosen k. 

3. Choose an initial non-singular, hermitian matrix, W 8 . Perform the numerical 
integration to compute the T-operator in (12.101) . 

4. Set the new h a ? to be h"? = (T^)" 1 . 

5. Return to item 3 and repeat until h a/3 approaches its fixed point. In practice, 
this convergence occurs in less than 10 iterations and does not depend on the 
initial choice of h a ". 

We now turn our attention to how one can measure convergence in Donaldson's algo- 
rithm. That is, for each finite step k, how does one define the error measure which 
tells us how close we are to the unique Ricci-flat metric? 



2.2 Ricci Curvature and Error Measures 

We begin by defining the Calabi-Yau volume form and several other useful quantities. 
Let X be a smooth Calabi-Yau variety of dimension d, and let Q denote the unique 
(up to a constant) (d, 0) form. The associated volume is given by 

Vol C y = QATi . (2.12) 
Jx 

Given an ample line bundle C, let denote the Kahler form 

= -j9§^i A % (2-13) 
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of the balanced metric associated with C k . Note that the Kahler class associated with 
this Kahler form is [uk] = 2nci{C k ). The associated volume is 

Vo\ k = j l J x UJ d k , (2.14) 

where u d denotes the (d, d) volume form u A • • • A u. As expected, when numerically 
computed, this value is close to k d J x uf. If, for example, we restrict to the case of 
the degree d + 2 hypersurface in P d+1 , then by Poincare duality we have that 



x «,*=(d + 2) I (uT +1 Y, (2.15) 



which allows us to fix the overall normalization to j x uf = d + 2. 

These two volumes allow one to define several measures of the convergence of the 
balanced metric to the Ricci-flat metric. The first is a, introduced in [1BJ H9] and 
given by 



1 



n Afi/Vol C y 



dVol cy . (2.16) 



Vol C y Jx 

Note that the Kahler form in ( 12 . 1 3[) is the Calabi-Yau Kahler form if and only if 
uif. is proportional to A Q. Since we know the exact volume form Q A Q, only 
Uk is approximate in (12. 161) and depends on the degree of k. Hence, the integral in 
(12. 16ft vanishes if and only if Uk is the Calabi-Yau Kahler form. As k is increased in 
Donaldson's algorithm, should approach zero. In [4*8l |4*9] . there is a prediction as 
to how Ok should converge to zero. They find that the error approaches zero at least 
as rapidly as 

for some constants a^. This prediction is verified by our results in the next subsection. 

The second error measure we use is a global measure of convergence for the Ricci 
scalar, 

||i% = — / \R k \ dVol fe . (2.18) 

VolcY J 

Rk is the Ricci scalar computed with the balanced metric and we integrate its 
absolute value. Note that by including the factor Voljj/ d , the naive linear k dependence 
of the integral inherited from the fact that ci(C k ) = kci(C) is canceled. As the 
balanced metrics converge to the Ricci-flat metric for increasing k, this measure should 
approach zero. 

Third, instead of using the exact Calabi-Yau volume form, one can also use the 
volume- form computed from the approximate Calabi-Yau metric. This is nothing but 
the Einstein-Hilbert action, so we define 

EH k = Volt d)/d J Rk ^ = Volt d)/d J R k V^Tk d 2d x . (2.19) 
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Note, however, that on a Kahler manifold the closed two-form Ric = Rfj dz l dzP defines 
the first Chern class of X, since C\(X) = [f^f]- Furthermore, Ric Aw^" 1 = 2dRu d and 
hence, the Einstein-Hilbert action 

EH k = Vo\ { ^ d]/d [ R k u d = Volj; 1 -^ / ci(X) A = (2.20) 
Jx Jx dl 

vanishes on a Calabi-Yau manifold for any metric (and hence for any integer k). As 
a result, we can use this to directly test the accuracy of our numerical integration. 
However, a better measure of the convergence of the balanced metrics to the Ricci-flat 
metric is given by 

\\EH\\ k = Vol^^y \R k \ y^t^ k d 2d x . (2.21) 

On a Calabi-Yau manifold, ||-R||fc, H-E'-^'llfe — 0(k~ x ) as k — > oo and, hence, these error 
measures should approach zero. 



2.3 The Quartic K3 

In the next two subsections, we use the algorithm outlined above to numerically 
compute the Ricci-flat metric on two Calabi-Yau manifolds. Two separate computer 
implementations of Donaldson's algorithm were developed in |l8lB9] and [501 S3] . We 
have utilized these two independent sets of code to check the accuracy of our results. 
These will be denoted as Codel and Code2 respectively in the following. 

Begin with the simple case of the K3 twofold. In this subsection, we consider the 
one parameter deformation X of the Fermat quartic in P 3 given by 

3 3 

X: J2 z i ~^Y[zi = 0. (2.22) 

i=0 i=0 

Here Zq, . . . are the homogeneous coordinates on P 3 . The deformation parameter i/j 
is, on the face of it, a complex number. However, redefining 

^ i—). e 27ri />, [zq :z x :z 2 : z 3 ] H- [h- e " 2 ™/ 4 z : z x : z 2 : z 3 ] (2.23) 

leaves X invariant, leading to a Z 4 isometry acting on the naive modulus ip. In 
other words, the actual modulus is ip A , and the complex structure moduli space has 
an orbifold singularity at the Fermat point ip = 0. As is well known, the algebraic 
variety X is smooth as long as ip is not a fourth root of unity and is away from the 
large complex structure limit, that is, ^ 4 ^ 1, oo. 

The very ample line bundle we consider is £ = Ox(l)- For each k, we computed 
the quantities in question for both ip = 0.5 and ip = 0.5e 27 ™/ 4 . These two parameter 



12 



crfc for if; — 2, Codel + 
(Tfc for ^ = |i, Codel 
(7fc for ^ = |, Code2 X 
(Tfc for ^ = |i, Code2 
Fit for k > 3: a k = 0.90/ k 2 - 0.25/A; 3 
a k = l.0b/k 2 - 1.52/A; 3 + 2.52/A; 4 
o-fc = 0.31exp(-0.42A;) 



k = 



k = 2 



k = 4 



k = 6 



k = S 



k = 10 



Figure 1: The error measures a k defined in Subsection \2.%\ The data shown 
is for the Quartic K3 defined as a hypersurface in P 3 (|2.22p . The 
complex structure parameter is chosen to be if) = \ and if) = |. 
Shown is data generated using the code developed in \4ty \4^f (Codel ) 
and data generated by the implementation in jl50\, \43)f (Code.2). The 
error measure is fitted to the theoretical convergence given in (|2.17p . 



\\R\\k for i> = f , Codel + 
\\R\\k for ip = \i, Codel 
\\R\\k for "0 = 1) Code2 X 
\\R\\k for if) = \i, Code2 
Fit for k > 3: ||i?|| fe = 1.55/fc + 52.89/A; 2 
;, \\R\\ k = 3.25/Jfc + 36.32/fe 2 + 35.41/A; 3 



h — -ga 



k = k = 2 k = 4 k = 6 k = 8 k=10 

Figure 2: Average scalar curvature for the metric on K3 C 



b3 
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values are a significant check of the validity of our numerical implementation since 
these two K3 metrics are, in fact, identified under the Z4 discrete isometry of X. 
We would expect the metric results for these two values of ip to be identical up to 
numerical errors. This is confirmed by the data shown in Figures [H [21 and [31 

In these plots, the T-map was iterated with 1,600,000 points, and the error mea- 
sures were computed with 500,000 points. The high accuracy allows one to conclude 
that both code bases give the same result. Moreover, especially for high k the as- 
sumption that Ok = O(t) fits our data better than exponential fall-off, see [US, l4"6] . 



2.4 The Quintic Threefold 

Now consider a Calabi-Yau threefold; specifically, the one parameter deformation X 
of the Fermat quintic in P 4 given by 

4 4 

X : 4 ~ 5 ^ II Zi = • ( 2 ' 24 ) 

i=0 i=0 

The algebraic variety X is smooth as long as ip is not a fifth root of unity. The T- 
map was iterated with 2,000,000 points, and the error measures were computed with 
500,000 points. 

The very ample line bundle we consider is again C = Ox(l)- For each k, we com- 
puted the quantities in question for both tp = 0.5 and ip — 0.5e 27 ™/ 4 . As described in 
the previous subsection, due to the discrete Z5 isometry of the quintic given in (12.241) . 
we would expect the metric results to be identical for these two choices. Once again, 



this is verified by the plots in Figure 4 and [51 We also repeated the computation of 



compared to ||-R||fc- The result is shown in |Figure 6[ As discussed in the Sub- 
section [2721 using the identity in f)2.20p one can estimate the accuracy of the numerical 
integration from EH^ « 0. 



3 Hermitian Yang-Mills Connections 
3.1 Generalizing Donaldson's Algorithm 

As we saw in the previous section, Donaldson's algorithm is a powerful tool for numer- 
ically approximating the Calabi-Yau metric. In this section, we investigate a general- 
ization of these techniques which can be used to approximate the field strength, F^- 1,l \ 
of a holomorphic connection which satisfies (11.131) . As discussed in ISection 21 Don- 
aldson's algorithm can be viewed as a method for numerically obtaining a particular 
Hermitian structure on the ample line bundle C k , see (I2.5p . This "balanced" fiber met- 
ric on C k allows one to define a balanced embedding of the Calabi-Yau space X into 
pnfc-i_ By mapping the coordinates x G X into the global sections s a G H°(X,£ k ), 
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+ 



R\k + 



\R\\ k = 1.55/fc + 52.89/fc 



2 



ML + 



|M|| fe = 2.63/k + U.03/k 



2 



-* X * X * X *- 



A; = k = 2 k = 4 k = 6 k = 8 k=10 

Figure 3: Curvature measures for K3 C P 3 . These provide the additional 
checks of numeric accuracy described in \Subsection jjl 



a k for ip 
a k for ^ 

cr fe for ijj 
o~k for V 
Fit for k > 3: cr fc = 3.51/A; 2 - 
% <j k = 3.51//c 2 - 5.12/fc 3 + 



|, Codel 
±e 27ri/5 , Codel 
= |, Code2 

. „27Ti/5 

2 ' 



Code2 

5.19/A; 3 
-0.14/& 4 



+ 

□ 
X 

o 



A; = fc = l A; = 2 A; = 3 A; = 4 A; = 5 k = 6 k = 7 k = 8 k = 9 

Figure 4 •' The a k error measures for the quintic threefold ()2.24p in P 4 . Shown 
is the error measure described in Subsection \2.2[ evaluated for the 
two values if) = ^ and ip = |. Codel and Code2 are associated to 
the implementations of fJSj and \UIk \4~3Tj respectively. 
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\\R\\ k for ip = i Codel + 
||i?|| fc for ijj = |e 2 ™/ 5 , Codel 

||i2||fc for ij) = i, Code2 X 
H^llfc for V = |e 2 ^/ 5 , Code2 
HtVor fc > 3: ||U|| fc = 18.84/fc + 71.29/A; 2 
S ||i%. = 0.13//C + 247.15/A; 2 - 368.23/A: 3 



fc = fc = l A; = 2 A; = 3 fc = 4 fc = 5 & = 6 & = 7 A; = 8 k = 9 
Figure 5: Average scalar curvature for the metric on the Quintic, (|2.24|) . 



+ 
+ 



+ 
+ 



\\R\U + 

\\R\\ k = 21.04/A; + 63.67/A; 2 

\\EH\\ k + 

\\EH\\ k = 20.20/*; + 52.00/A; 2 

EH k X 



X * X X X X X 





k = k = l k = 2 k = 3 k = A k = 5 k = 6 k = 7 k = 8 
Figure 6: Curvature measures for the Quintic, (|2,24D . 
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that is, 



x h 



[s (x) 



(*)] 



(3-1) 



we produce a map i\~ : X — » P nfc_1 = G(l,rife — 1) where nk = h°(X,C k ). The pull- 
back of the associated Fubini-Study metric was shown in ISubsection 2.11 to converge 
to the Ricci-flat metric on X in the limit that k — > oo. It is a natural question to ask 
whether or not an analogous algorithm could be developed to approximate Hermitian 
metrics on higher rank vector bundles? In particular, could one find an approximation 
scheme to produce an Hermitian metric on an arbitrary stable bundle V of rank n 
such that it satisfies condition ( 11. 81) ? Fortunately, precisely this question has been 
addressed in the mathematics literature [12] and in [35] . 

To begin a generalization of Donaldson's algorithm, consider defining an embed- 
ding via the global sections of a twist of some holomorphic vector bundle, V, with 
non-Abelian structure group. That is, consider a map 



x h 



fSl{x)\ 
\S?(x)) 



( S k( x )\ 



(3.2) 



from the coordinates, x, of X into the global sections S* £ H°(X,V (8) C k ), where 
a = . . . N k — 1 = h°(X,V (g) C k ) (the number of global sections of V <g> C k ) and 
the index a — 1, . . .n is valued in the fundamental representation of structure group 
K C ?7(n) of the rank n bundle V. We hope then to define the embedding 



X -»• G(n, iV fc - 1) 



(3.3) 



where G(n, Nk — 1) denotes the Grassmannian of the relevant dimension. By the 
Kodaira embedding theorem [61], given a holomorphic vector bundle V and an ample 
line bundle C, there must exist a finite integer /c such that, for any k > k , a "twist" 
of the bundle V(k) = V®C k defines an embedding, : X — > G(n, N/. — 1). 

Such twisting is a clear necessary for the bundles of interest in E% x Eg heterotic 
compactifications. Because a slope-stable bundle V in a heterotic compactification 
cannot have any global section^, H°(X,V) = 0, we begin by twisting the bundle by 
some sufficiently large power of an ample line bundle C. That is, we consider V ® L k . 
If C is ample, then V <g> C k will be generated by its global sections; that is, it will 
define an embedding as in ( 13. 31) . In our search for a solution to the Hermitian Yang- 
Mills equation ( I1.13p . the connection on the twisted bundle will be closely related 
to the original connection, since such a twist only modifies the trace part of the field 



1 To see why this is the case, note that for the structure group K of V to embed inside Eg one 
must require that K = SU(n) for n = 3,4, 5 and, hence, that C\(V) = 0. Therefore, from (|1.10[) it 
follows that /i(V) = 0. However, if H°(X, V) ^ then Ox must inject into V in contradiction with 
the assumption of stability. 
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strength. Stated in terms of algebraic geometry, the process of twisting will not modify 
the slope-stability properties of V since V <8> C k is stable if and only if V is. 

As at the beginning of ISubsection 2.1] where we chose the trial form of the Kahler 
potential in (j2.4p . here we begin with another simple anzatz; this time, however, for 
the Hermitian structure G in (jl.4p . We consider an Hermitian matrix (G~ 1 ) ab of the 
form 

N k -1 

(G- X ) al = J2 Ha ^ a a {S% , (3.4) 

where H a/3 is an arbitrary matrix and are the global sections of V <g> C k . As in 
( II. 9p . this fiber metric induces an inner product on the space of sections, S^, given by 

N f - N f ___ __ 

<Sf,\S a >=— / S^G^-'SUVoIcy = / S^H^Slr'SUVolcv . 

Voicy Jx Volcy Jx 

(3.5) 

With this definition of the inner product on sections, one can give a natural gen- 
eralization of the T-operator (12.10p . This generalization, 

T(H) a0 = J S a (&HS) ^Sp dVol cy - , (3.6) 



was introduced in (42] and studied numerically in (43]. Note that when V is a line 
bundle, (13.61) reduces to (12. 101) and we return to the case of a balanced embedding into 
P^- 1 . As in the previous section, we will describe how the iteration of the generalized 
T-operator can produce a fixed point which describes an Hermitian-Einstein bundle 
metric. 

In [12], it was shown that the bundle V is Gieseker stabl __ if and only if the fc-th 
embedding (defined by V ® C k as in (13. 3p ) can be moved to a "balanced" place. That 
is, if there exists an orthonormal section-wise metric on the twisted bundle such that 

(T{H) aB y l = H** (3.8) 

is a fixed point of the generalized T-operator We can use this special metric on V ® C k 
to define an Hermitian metric on V itself. Let Gc denote the Hermitian metric on C 
and G^ the balanced metric on V ® C k . Then 

G k = G {k) ® G~ c k (3.9) 

is an Hermitian metric on V. This appears in the following important theorem [4~2] . 



2 Let C be an ample line bundle. For any torsion-free sheaf F define the Hilbert polynomial with 
respect to C as 

Pc(F){n) = 7-—T- (3.7 

rank{J- ) 

where ^(J 7 <E) C n ) is the index of T <S) C n . Given two polynomials / and g, we will write / -< g if 
f(n) < g(n) for all n 3> 0. Then a bundle V is said to be Gieseker stable if for every non-zero torsion 
free subsheaf JcV, p C (T) -< p £ (V) [54]. 
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Theorem 4 (Wang). Suppose V is a rank n, Gieseker stable bundle. If Gk —> Goo as 
k — > oo, then the metric G^ solves the "weak Hermite- Einstein equation" 

9 fj F fj =(fi+ ^y^) Uxn (3.10) 

where 

• R is the scalar curvature. 

• R = J Ry/ det g d 2d x. This vanishes (for any Kahler metric) on a manifold with 
vanishing first Chern class. 

Procedurally, the process of obtaining the Hermitian-Einstein fiber metric on a 
slope-stable bundle, V, is very similar to that outlined for the Ricci-flat connection 
in Section [2] - for each value k of the twisting, we iterate the T-operator associated 
with the embedding defined by H°(X, V x C k ) until a fixed point is reached. Then, by 
Theorem HI in the limit that k — > oo the induced connection solves (I3.10p . However, 
there is an immediate and important difference between this generalized algorithm 
and Donaldson's algorithm for Ricci-flat metrics. Note that while all slope-stable 
bundles are Gieseker stable [51] . the converse does not hold-not all Gieseker bundles 
are slope-stable. That is, there certainly exist cases where the iteration of the T- 
operater does not converge (at a given k). However, if V is a slope-stable holomorphic 
bundle, then the iteration H n+ i = T(H n )~ l will converge at each k, and in the limit 
that k — > oo, produce the Hermitian bundle metric satisfying 113. 10p (via its 
associated field strength defined in (11.71) and fll.Sp ). Moreover, in the case where the 
Calabi-Yau metric g^ is Ricci-flat, (13.101) simply reduces to fll . 13j) . Thus, we have 
found a solution to the standard Hermitian-Einstein equations. 

3.2 Untwisting 

Despite having described a balanced embedding associated with the twisted bundle 
V <g) C h , and an Hermitian metric G^ satisfying 1 13. 10p . our task is not yet complete. 
Thus far, the discussion has been completely general and could be applied to any 
stable U(n) bundle. However, for the physics associated with an E$ x Eg heterotic 
string theory, we are ultimately interested in the explicit connection, A, on the SU(n) 
bundle V itself satisfying the slope-zero Hermitian Yang-M ills equations, (11. ip . Since 
the process of twisting V by a line bundle C k in the above construction clearly modifies 
the trace-part of the connection, one must subtract this line bundle contribution to 
get the connection on V only. To do this, we have to separately compute the balanced 
metric, Gc, on L. 

Each of the bundle metrics, Gc on £ and G^ on V <g> C k , as well as the Calabi- 
Yau metric g 1 ^ must be approximated numerically. The computation of the Hermitian 
Yang-Mills connection in fll.ip rests on three finite-dimensional approximations: 
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1. The degree k g at which we compute the metric, 

2. The degree kn used to compute the twisted connection 

3. The degree kh used to compute the U(l) part of the connection that must be 
subtracted to obtain the final SU(n) connection on V. 

Since the Hermitian Yang-Mills connections are unique (for a given choice of Calabi- 
Yau metric and geometric moduli), if we can approximate the connections on V ® 
C ku and C kh with sufficient accuracy, then we can compute them numerically in a 
completely independent way. In terms of the metrics on V® C kn and C kh , we must 
find 

• the metric on V<g>£ fcff , G^ kH ' = (S^HS), where kn denotes the finite number of 
iterations that will be performed to compute the connection. 

• the metric on L kh , (Gc) kh = (s^hs) (and hence Gc = (s^hsY^ kh ), where kh will 
denote the number of iterations to approximate the connection. 

Then, as in (13 .9D we find that 

G = G {k « ] x G~ c kH = (S^HS) (s^hsY kH/kh (3.11) 

is the fiber metric (TOD on V. As before, S G H°(X, V <g> C k ") and s G H°(X, C kh ). 

Using (I1.7P and (II. 5J) , in terms of the Hermitian metric the connection on V is 
then simply 



A(V) = d 



( S^HS) ( s^h 



(ViLs) (s^h 



_ k H /k h 
S ' 



[d{S^HS]) (S*HS)~ L - ^(d(jhs)) (sW) ' (3.12) 



= A(V®L kH ) - ^A(L kh ) . 

k h 

That is, one can "untwist" the connection simply by subtracting the trace of the U(n) 
connection to produce an SU(n) connection. The degrees of approximation kn and 
kh can, in principle, be chosen independently, as long as both are sufficiently large. 
We will discuss the choice of the line bundle £ and twisting degree kh in more detail 
in the next subsection. 

The curvature is given by 

F (0,2) = F (2,0) = g fj Ff , = g i] d]Ai = g ijg ]di ln (tf HS \ ^ hs ^j ~ kH/k \ (3 13 ) 

In summary, the following is an outline of the computation the Hermitian Yang- 
Mills connection (that is, the solution to (II. ip ). 
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1. Following Donaldson's algorithm derived in ISubsection 2.11 approximate the 
unique Ricci-flat metric g l i on X to the desired degree (determined by k g , the 
iteration parameter of Theorem 2, (12.91) ). 

2. For a given holomorphic vector bundle V, choose an ample line bundle £ on X 
and a degree k# (that is, a twisting of the vector bundle V (g) C ka ) at which to 
compute the "balanced" fiber metric (13.8j) 

3. Calculate a basis {S'q,}^^ 1 for H°(X,V <8> £ fcff ) at the chosen 

4. Choose an initial, invertible, Hermitian matrix, if 7(5 for the ansatz (I3.4p . Per- 
form the numerical integration to compute the T-operator in (13. 6ft . 

5. Set the new H a P to be = (T^)" 1 . 

6. Return to item 3 and repeat until approaches its fixed point. In practice, 
this convergence occurs in less than 10 iterations and does not depend on the 
initial choice of H a ^. At this point, we have obtained the Hermitian metric on 
V <S> C ku to the desired accuracy. 

7. Approximate the Hermitian metric on £ to a chosen degree kh by repeating 
steps 2 to 6 for C kh 

8. Compute the "untwisted" connection and field strength via A3. 12[) and (I3.13P 
We turn now to a discussion of the accuracy of these numerical approximations. 



3.3 Measuring the Error 



As discussed above, in a computer implementation of the generalized Donaldson's 
algorithm one must rely on three expansion parameters. Two are associated with the 
approximation of the balanced Hermitian metrics on V® C kn and C kh , and the third 
with the Calabi-Yau metric calculated with parameter k g . Having implemented this, 
we would like to know how far our numerical connection for fixed k g , kn, kh deviates 
from the exact Hermitian Yang-Mills connection (II. ip . 

F 0,2 and F 2,0 are automatically zero with our ansatz (I1.7P for the connection. 
Hence, for an SU(n) bundle we have only to test how far the "color matrix" g^F^ 
deviates from the zero matrix. Note that the matrix entries depend on the chosen 
frame. Thus the appropriate invariant quantity to consider for such an Hermitian 
matrix is its (real) eigenvalues 



\ 



(3.14) 



A, 
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For an SU(n) bundle, we expect 



Aj V i (3.15) 

if and only if the connection is approaching a solution to (11. ip . We begin by investi- 
gating the eigenvalues of the "color matrix" at a point. An integrated error measure 
will be discussed in the next subsection. To illustrate these concepts, we turn to a 
simple example of a stable bundle. 



3.4 An Example 

To begin, consider the Quartic K3 of ISubsection 2.31 We define an explicit such 
hypersurface in P 3 by 

|4 + z \ + 4 + ^ _ 2*o W 3 = 0} C PjU i: ^ 3] , (3.16) 

where Z{ denote the coordinates of P 3 . The Ricci-flat metric on this manifold was 
computed up to degree k g = 10 in ISubsection 2.31 and the results plotted in Figures 
HJiandEJ 

On this Calabi-Yau twofold, we define the following rank 2, holomorphic vector 
bundle with structure group SU(2). This sample bundle is defined through the so- 
called "monad" construction [SH E21 ESI 125] 

— ► 0(-3) A O(-l) © 0(-l) © O(-l) — ► V stablc — ► . (3.17) 

Here V stable is defined as the cokernel of the map / = {xq,x\,x%) between the direct 
sums of line bundles. Using the techniques of [M| 1(32"] . it is straightforward to prove 
that V stable in (I3.17P is slope-stable. Hence, by ITheoreni~4t we expect the T-operator 
to converge. 

As discussed in ISubsection 3.2| to apply the generalized Donaldson algorithm one 
must define the embedding : X — > G(n,Nk H ). To do this, we must compute the 
global sections of the twisted line bundle V © C ku for some ample line bundle C. 
Fortunately, for the bundle defined in ( 13. 17ft the global sections H°(X, V © C kli ) can 
easily be computed as follows. Letting C ku = 0(m), from (I3.17P we induce the twisted 
short exact sequence 

— ► 0{m -3)-?-> 0(m - l)® 3 — > V stable © 0{m) — > . (3.18) 
Then the global sections are given simply as the cokernel 



22 




where both parts of this quotient are the global sections of sums of ample line bundles 
when m > 3. Furthermore, on the Quartic K3 the global sections of the line bundle 
0(n) can be computed as a simple polynomial space of dimension 

n < 

< n < 4 (3.20) 
"(V) "> 4 

With these definitions in hand, one can compute a basis of polynomials of H°(X, V eg) 
C kn ) of the degree given by (13.201) . 

Following the algorithm developed in the preceeding sections, we first compute 
g^Ffj at the randomly chosen point 

P = [1 : 0.707124 + 0.707124z : 0.1 : 0] . (3.21) 

The bundle V is rank 2, so there are two eigenvalues Ai,A2 in f)3.14p . Define the 
following extremal quantities. 



|A| max = max {|Ai|, |A 2 |} 
|A| min = min {|Ai|, |A 2 |} 



(3.22) 



As we iterate the twisting degrees ku and kh associated with the connections on 
V ® C k " and C kh respectively, we can plot the eigenvalues of (13.141) at the point 
( 13.211) . For this example, the connection T-operator was iterated 7 times, numerically 
integrating over 1,000,000 points. As expected, we find that A^ — > as kh, ku — > oo. 



The results are shown in |Figure 7| Inspecting this graph, we note several important 
features. First, observe that the eigenvalues will converge to zero along any diagonal 
ray, that is, in any limit as kh,kn —> oo. Second, note that there is clearly some 
variation in the rate at which the eigenvalues along these rays converge to zero. We 
will define the optimal choice in the next subsection. 



3.5 Subtracting the Trace 

As we saw in ISubsection 3.2"! after twisting up the vector bundle to V <S> £j ka there is 
a choice in how we compute the connection on the line bundle C As long as there 
is some limit (for example, kh — > oo and picking the balanced connection) where 
the connection on the line bundle approaches constant curvature, we will eventually 
approximate the slope-zero Hermitian Yang-Mills connection on V. Of course, some 
choices of line bundles work better than others. In this subsection, we discuss the 
optimal choice-the determinant line bundle of V <g> L kH with the induced connection. 
The determinant line bundle is defined as A n (V <g> L kH ). Taking this line bundle as the 
"untwisting" line bundle and choosing kh = i&nk(V)kH, we then choose the Hermitian 
metric 

(G c ) kh = det(G {kH) ) = det (S^HS) (3.23) 
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Figure 7: The eigenvalues of g tJ F i j(P). The actual solution to the Hermitian 
Yang-Mills connection is characterized by having A = at each point. 
We see that the numerical solution approximates this in the limit as 
ku, hi ->■ oo. 
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on (A n (V® L kH )) k K 

Let A (/ch) be the "twisted" eigenvalues of g fj F^ H) on V <g> £ fc ", and let A be the 

corresponding eigenvalues of g^Fq on V after "untwisting". Then, from (13 . 1 3[) we 
have ^ ^ 

A. t = \f H) tr F$ b) = \f H) - ■ (3.24) 

rankV y w 4 rankV v ; 

Therefore, the effect of this untwisting is precisely to subtract, at each point, the 
average of the eigenvalues. In this sense, using the determinant line bundle is the 
optimal way to untwist. Furthermore, with this choice of untwisting, the approximate 
connection on V is SU(n) at each value of kh (as opposed to the other choices for 
"untwisting" line bundles for which it is only as kh becomes large that the connection 
approaches SU(n), see Figure [8]). In Figure 8, we compare the untwisting done with 
various line bundles C kh = 0(1) h = 0{kh) on the Quartic with the untwisting by the 
determinant line bundle. The determinant line bundle clearly produces the most rapid 
convergence to zero eigenvalues. Therefore, for the remainder of this paper, regardless 
of the choice of "twisting" line bundle, C, we will always use the determinant line 
bundle A n (V <g> C kn ) to untwist. 



3.6 Integrated Error Measure 

To define a true error measure for the approximation to the Hermitian Yang-Mills 
connection, we must test the approximation at all points (as was done for the error 
measures for the Ricci-flat metric in lSubsection 2.21) and, hence, integrate over X. To 
do this, define the error measure 

t(Av) = ~. / (£WW M * . (3.25) 

2vr Voh 9 rank(V) J x \^ J 

This is simply a global check of the eigenvalues in (13.1 4ft . To understand the nor- 
malization in (13.251) . first note that scales as ^- with the degree k g of the metric 
computation. This explains the prefactor k g in ( 13.251) . We also divide by the metric 
volume Voh g to cancel the scaling of the volume element y/g. Finally, observe that 

• If all eigenvalues are positive, then r(Ay) = f Ci(V) A u D ~ l G Z. This is the 
reason for the prefactor since [u>k] = 2nci(C k ). 

• If C\(V) = 0, then Tk H (Ay) — > measures convergence to the Hermitian Yang- 
Mills connection. The eigenvalues must occur with both signs. 

Since Ci(V <g> C k ») = c x (V) + nk H ci(C) is the first Chern class of V <g> C k " (with 
n = rank(V)), we can predict the values of rj: w ^ for the twisted bundle V eg) C kn as 
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Figure 8: The eigenvalues of g %3 F^{P) of the rank 2 bundle V on the Quar- 
tic for different choices of "untwisting line bundle" C kh = 0(l) kh 
(same data as in Figure 7\ ) as well as the optimal untwisting with 
the determinant line bundle. The upper and lower boundaries of the 
colored bars are \X\ max and |A| m j ra , respectively. In the determinant 
line bundle case, the two eigenvalues are ±A and their magnitude is 
exactly equal at each value of kh- This magnitude is shown by the 
blue crosses. 
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we increase kjj. Hence, we can learn something from evaluating the r-integral for the 
twisted as well as for the untwisted connection. Define 



(tw) 



A'// 



~{A V(g)L k H ) 



(3.26) 



where l ra nkv is the identity matrix in the gauge indices. That is, rj^ is the (properly- 
normalized) integral over the eigenvalues of the twisted-up bundle and t^ h is the 
integral over the eigenvalues after untwisting with the determinant line bundle (this 
expression follows immediately from (I3.12p by taking C kh = (A n (V <g) £ k H^yank{V)k H y 
Note that it implicitly depends on the degree k g at which we computed the metric. 
In ITable 1\ we list both and r kiI for the rank 2 bundle defined by (13.171) in 



k H 


/ Ci(V® C kH ) Alo 


(tw) 




2 


4 


4.00057 


0.904183 


3 


6 


6.00329 


0.765671 


4 


8 


7.99855 


0.532639 


5 


10 


9.99944 


0.331912 


6 


12 


12.0012 


0.229044 


7 


14 


13.9983 


0.170002 


8 


16 


15.9989 


0.132257 



Table 1: The integrated error measure for the stable bundle (|3.17j) from 



Subsection 3.4- T [ before untwisting computes the integral first 



Chern class. The comparison between the exact value for the inte- 
gral Chern class and t£™ is made in the first two columns. This 
comparison provides a confirmation of the validity of the algorithm, 
even before untwisting. In the final column, the error measure for the 
untwisted bundle is presented. 

ISubsection 3.41 In the first column, the integral over the first Chern class of V ® L kH 
is listed, where C = 0(1) on the K3. As expected, rf^' closely matches this value. 
In the last column, the untwisted error measure is presented. It clearly approaches 
zero as kn increases. 



In |Figure 9| we present our first explicit example of the convergence of the gener- 
alized Donaldson algorithm for a slope-stable bundle. We have explicitly computed 
the field-strength of the SU(2) bundle in (I3.17p . We find that the associated inte- 
grated error measure, Tk H in (13. 25 p . is converging to zero in the limit that kn » 1. 
In addition, in Figure 9] the form of the decreasing t^ h values is compared with the 
prediction of an 



1/kff fall-off in the error measure predicted in 
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Figure 9: Tk H for the stable bundle V 'stable defined in (|3.17j) on the Quartic K3. 

4 Stable vs. Unstable 

In this section, we investigate the behavior of the generalized Donaldson algorithm 
applied to arbitrary SU (n) vector bundles. This behavior will exhibit more complexity 
than we encountered in the computation of the Ricci-fiat metrics of ISection 2l In 
applying the Donaldson algorithm to the computation of the metric, we are guaranteed 
that for any compact Kahler manifold X with c±(X) = the Ricci-fiat metric exists. 
Furthermore, the algorithm presented in lSubsection 2.ll will approximate it in the limit 
that kg — >■ oo. However, the situation is different for vector bundles. A Hermitian- 
Einstein metric does not exist for every vector bundle V on X with Ci(V) = 0. As 
discussed in Sections [T] and 13.11 the vector bundle must be slope poly-stable for such 
a Hermitian-Einstein metric to exist. Since this property is difficult to guarantee 
initially, it seems prudent to ask what behavior we expect to see if we apply this 
algorithm to bundles that do not, in fact, admit an Hermitian-Einstein fiber metric? 
In the following, we demonstrate that, surprisingly, the Donaldson algorithm produces 
distinctive and interesting results even in the case where the bundle is not stable. 
Furthermore, we will argue that in view of the difficulty in proving slope-stability 
of vector bundles, the generalized Donaldson algorithm provides a new and attractive 
way of numerically deciding whether or not a bundle is stable. 

With this goal in mind, in the next sections we will investigate the slope-stability 
properties of vector bundles and how this behavior appears in the results of the gen- 
eralized Donaldson algorithm presented in ISubsection 3.11 
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4.1 Taxonomy of Slope-Stability 



Recall that the slope of a vector bundle V is defined, with respect to a given Kahler 
form to, as 



As discussed in ISubsection 1.31 there are four possible types of behavior for holomor- 
phic vector bundles. These are [51] 

1. Stable: An indecomposable bundle, V, is called stable if, for all subsheaves 
T C V, /i(JF) < fx(V). 

2. Poly-stable: A bundle is called poly-stable if V is a direct sum of stable bundles 
with the same slope: V = i Vj with /x(Vj) = /x(V) V z. 

3. Semi-stable: A bundle is called semi-stable if for all subsheaves J 7 C V, /i(J-") < 
fi(V). It is properly (or strictly) semi-stable if V is semi-stable but not poly-stable 
(that is, it is indecomposable). 

4. Unstable: An unstable bundle is one for which ^(J 7 ) > /x(V) for at least one 
proper subsheaf J 7 C V. Unstable bundles can be either decomposable or inde- 
composable. 

Recall that a Hermitian-Einstein metric exists only for the first two entries in this 
list, that is, for poly-stable bundles (stability being a subset of poly-stability). 

Before we discuss the results of the generalized Donaldson's algorithm in the four 
cases above, we must make one further observation about slope-stability. The pre- 
ceding definitions provide a measure of the degree of substructure present in a given 
bundle. More precisely, stable bundles are simple |54j . that is, H°(X,V (g) V v ) = C 
and any morphism between two stable bundles Vi,V 2 with rank(Vi) = rank(V 2 ) and 
Ci(Vi) = c 1 CV 2 ) is an isomorphism. Likewise, poly-stable bundles are direct sums 
of simple objects. But properly semi-stable and unstable objects exhibit a richer 
structure. Theoretically, there is a natural framework for describing the behavior 
of semi-stable and unstable vector bundles in terms of a decomposition into simple 
objects. We discuss this in the next subsection. 

4.2 The Harder-Narasimhan Filtration 

The following theorem makes explicit how unstable sheaves may be described in terms 
of semi-stable sheaves, and semi-stable sheaves in terms of stable sheaves |54j . 

Theorem 5 (Harder-Narasimhan). Given a holomorphic bundle V over a closed 
Kahler manifold X (with Kahler form u), there is a filtration (called the Harder- 
Narasimhan filtration) by subsheaves 




(4.1) 



= Jo C jF x C . . . Ti 



V 



(4.2) 
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such that TijTi-x are semi-stable sheaves for i — 1, . . . I and the slope of the quotients 
are ordered 

/iOFx) > fifo/Fi) > ■ (4-3) 

If V is semi-stable, then there is a filtration by subsheaves (called the Jordan-Holder 
filtration), 

Q = T Q <zFx<Z...F l = V (4.4) 

such that the quotients TijTi-x are all stable sheaves and have slope \i{Tii 'J-'i-x) = 
/i(V). In addition 

Gr(V) = T 1 ® T^jTx © • • • Ti/Fi-i (4.5) 
is uniquely determined up to isomorphism (and is called the graded sumj. 

One consequence of lTheorem 51 is a description of the moduli space of semi-stable 
sheaves. Two semi-stable bundles V\ and V 2 are called S-equivalent if Gr(Vx) = 
Gr(V 2 )- The concept of S-equivalence arises when trying to define the notion of a 
moduli space of sheaves. Stable bundles correspond to unique points in their moduli 
space, while a moduli space of semi-stable sheaves can only be made Hausdorff if 
each point corresponds to an S-equivalence class. It is worth observing that each 
S-equivalence class contains a unique poly-stable representative (namely the graded 
sum (1431) ). 

With these classifying notions in mind, we now return to our discussion of the 
generalized Donaldson algorithm and the results of the numerical scheme outlined in 
ISubsection 3.1[ If we numerically approximate the color matrix g^Fq for the four 
types of bundles described in points l)-4) above, what will we find? What will be the 
eigenvalues in (13.141) described in ISubsection 3.3| or the behavior of the error measure 
Tk H defined by (13.251) in ISubsection 3.61 ? 

To answer these questions, we again consider the four types of bundles described 
in ISubsection 4.11 above. Depending on the slope-stability properties of the bundle, 
the Harder- Narasimhan filtration indicates that there exists a connection which will 
produce a color matrix g %i ' Fq with the following behavior. 

1. If V is a poly-stable (including properly stable) bundle with slope /i, then 

/A* \ 



MV)l nxn • (4.6) 



By lTheorem 4\ in this case we expect that the generalized Donaldson algorithm 
will approximate a solution to the Hermitian Yang-Mills equation (jl.2p . 

2. If V is semi-stable, then, in general, it is not a solution to the Hermitian Yang- 
Mills equations. The form of its color matrix will depend on the values of the 
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bundle moduli, ^{EndiV)). These moduli determine where V is chosen to be 
within it S-equivalence class. In general, the color matrix g^F^ will produce 
constant, non-equal eigenvalues. However, by varying the bundle moduli, V can 
be made arbitrarily close to the unique poly-stable representative (the graded 
sum (14. 5p ) in its class. That is, for a properly semi-stable SU(n) bundle, de- 
pending on the choice of bundle moduli, the eigenvalues can be made arbitrarily 
close to zero. 

3. If V is unstable, there are two possible ways in which the field strength fails to 
satisfy the Hermitian Yang-Mills equations. Either 



(a) Its Harder-Narasimhan filtration f!4.2[) consists of a sum of bundles (locally 
free sheaves with constant rank) with inequivalent slopes, In this case, 
the color matrix takes the form 

1^2 



V 



(4.7) 



(b) Or its Harder-Narasimhan filtration ( 14. 2 p contains torsion-free sheaves (whose 
rank can change over higher co-dimensional loci in the base X). In this 
case 

/ fJ>i{x) \ 



(4.8) 



where fii(x) can vary in magnitude over the base X. In particular, at the 
locus in the base where the rank of V jumps, Fn can diverge to produce a 
curvature singularity. 

How will this behavior appear in the generalized Donaldson algorithm? As stated 
above, we are guaranteed that the algorithm produces the physical connection only 
in the case when the vector bundle is stable and the T-operator converges. But we 
can now ask, what happens when we apply the algorithm to properly semi-stable 
or unstable bundles? We will explore this experimentally first by looking at sample 
SU(n) bundles defined on the familiar Quartic K3. 



4.3 Examples of Semi-Stable and Unstable Bundles on the 
Quartic K3 

We have previously seen in lSubsection 3.4l the behavior of a slope-stable bundle under 
the generalized algorithm. For the stable bundle defined in (J3.17]) . the algorithm 
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produces convergence to an Hermitian Yang-Mills solution of the form ( 14. 6ft . The 
integrated error measure Tk H was shown to converge to zero in Figure 9 We now 



perform a simple experiment to see what happens when we apply the generalized 
Donaldson algorithm in a case where ITheorem 41 does not apply, that is, for a non- 
stable bundle. 



4.3.1 A Semi-Stable Bundle 

First, let us explore the behavior of a properly semi-stable SU(n) bundle. Such 
a bundle is not a solution to the Hermitian Yang-Mills equations, but can come 
arbitrarily close to a solution as we vary the bundle moduli. To illustrate this, consider 
a rank 3 bundle over the same Quartic K3, ( I3.16p . We choose the SU (3) monad bundle 

, (_4) Acm - Stablc : 0(-2) © O(-l)® 2 © O V scm , stablc , (4.9) 



where 



/ 3uz+2xz+7yz \ 

/'scmi-stablc 



(6u 2 x+20u 2 y+15uxy+5x 2 y+3uy 2 + 10xy 2 + 18u 2 z+21uxz 

+ 7x 2 z + 6uyz + 7xyz + 7y 2 z+8uz 2 + 16xz 2 + 10z 3 ) 



(4.10) 



i 7u 2 x+10uxy+4 : y 3 +6uxz+5uyz+7uz 2 +2xz 2 +10yz 2 +z :i i 

\ 8u 4 +6u 2 yz+5uy 2 z+5uyz 2 +19xyz 2 +7uz 3 



and x,y,z,u are the homogeneous coordinates on P 3 . Here the map, /scmi-stablc is an 
element of H°(X, Hom(A, B)), where A = 0(-4) and B = O © 0{-2) © O(-l)® 2 . 
The global sections of V semi _ st abie can be computed as in (13. 19[) and (13.201) . 
The Harder-Narasimhan filtration of this bundle takes the simple form 

0@F, (4.11) 

where J 7 is a rank 2 bundle with Ci(F) = defined by 

— >■ 0(-4) -)• 0{-2) © C(-l)® 2 . (4.12) 

This is an SU(2) bundle, but it cannot solve the slope- zero Hermitian Yang-Mills 
equations. Instead, as described in point 2 above, for a generic value of the bundle 
moduli (I4.10p we expect t^ h to converge to a constant, finite, positive number. 



4.3.2 An Unstable Bundle with a Filtration by Line Bundles 

As our next case, consider an unstable vector bundle with a Harder-Narasimhan 
filtration that is a sum of line bundles. We select a simple example for which V unsta bic 
itself is the following sum of line bundles over the K3 in (I3.16p . 

^unstable sum = 0(1) © 0(-l) . (4.13) 
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Here V is manifestly slope-unstable since /i(C(l)) > fi(y unsta bie sum)- From point 3(a) 
above, we expect the color matrix to take the form 




(4.14) 



This form is the same globally, since the bundle is a direct sum. As a result, because 
error measure T\. H in (I3.25P contains the sum of the absolute values of the eigenvalues 
of the color matrix, we expect that as kn — > oo error measure Tk H ~ |Ai| + | A2 1 — > 2. 



4.3.3 An Unstable Bundle with Sheaf Filtration 

Next, consider the case of an unstable bundle with a more complicated Harder- 
Narasimhan filtration. On the Quartic K3 given by (I3.16p . we define the unstable 
SU(2) monad bundle V unstable by 

— ► 0(-3) -A Oil) © G(-2)® 2 — ► V unatable — ► , (4.15) 

where 

/ = (x\x2+5xix\+10x'lx\+7x2x\+<ixQx\ ) Xi + 15x2+6x3, 2xo+x 2 +3xi + 12x3) (4.16) 

A simple analysis along the lines of [52] (see, for example, Hoppe's Criterion), reveals 
that V ung t a bie is a slope-unstable bundle. It is de-stabilized by a rank 1 sheaf, T . 
Furthermore, the Harder- Narasimhan filtration of V unsta bi e is given simply by 

Ox(-l)©J, (4.17) 

where J 7 is a rank one sheaf described by 

— > J 7 — > C(2)® 2 0(3) — > (4.18) 

with J 7 = ker(g) and ci(J-") = 1. This rank 1 object is a sheaf and not a line bundle 
because its rank can jump (from rank 1 to rank 2) over a higher co- dimensional locus in 
the base X. This occurs when the defining polynomial map, g = {x\ + . . . , 2xo + . . .) 
shares common zeros with the defining polynomial of the K3 surface in (12.22]) . It 
should be noted that V unstable itself has no such singularities. It is a vector bundle, 
despite the presence of subsheaves such as T in (I4.18|) . As a result, we expect the 
connection to produce a color matrix of the type shown in ( 14. 8P ; that is, 

9 % ~ • (4.19) 

Listed above are possible expected forms for the color matrix of several non-stable 
bundles. We now ask: what will the Donaldson algorithm produce when applied to 
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these examples? What behaviour do we expect for the error measure Tfe H ? How will 
the presence of filtration-sheaves, such as (14.171) . manifest itself? Before we investigate 
these questions, it should be noted that the Harder- Narasimhan filtration of V unsta bi e 
is almost the same as the form of the previous example in (14.131) . That is, the rank 
1 sheaf T nearly everywhere resembles the line bundle 0(1) except at points. As 
mentioned above, had this been the case there would be a clear prediction for the 
results of the integrated error measure. But what happens for V unsta bic as defined in 
(14.151) ? We turn now to an examination of the color matrix eigenvalues. 

We present the numerical results for the three bundles above, as well as for the 
stable bundle (13.171) from ISubsection 3.41 The error measure Tk H is plotted for these 



bundles in |Figure W\ Despite the fact that these bundles are not solutions to the 
Hermitian Yang-Mills equations, we find that the way in which they fail to provide a 
solution are in exact agreement with the mathematical structure of semi-stable and 
unstable bundles. As expected from the Harder-Narasimhan nitrations, the three 



cases shown in Figure 10 can be distinguished by the behavior 



{0 (stable) 
const. > (semi-stable) (4.20) 

const. > or oo (unstable) 

The error measure associated wth the semi-stable bundle V sem i_ sta bi e converges to a 
non-zero constant value. Furthermore, the unstable sum of line bundles produces a 
constant value Tk H = 2. Somewhat surprisingly, we find that the error measure for 
^unstable converges to a constant value of ~ 2, as if its Harder-Narasimhan filtration had 
been 0(— 1) ©C?(+l) (that is, of the /i = const, type described in (14.71) ) instead of the 
actual filtration in ( 14.17p . To understand the result for this last unstable bundle, one 
must look in more detail at the computation of the field strength and the behavior of 
the subsheaf J 7 in the algorithm approximating the connection. This will be explored 
in detail in the next subsection. We will give an explicit discussion of how the T- 
operator fails to converge in this case, and how the generalized Donaldson algorithm 
none-the-less produces a connection that exhibits the correct physical singularities. 



4.4 Eigenvalues Along a Geodesic 

In this subsection, we explore in detail the behavior of the unstable bundle V unsta bi 
on the Quartic K3. As discussed above, while the rank of V unsta bi e remains constant 
everywhere on the K3 surface X, the rank of a de-stabilizing subsheaf J 7 in (14. 18ft 
can increase over a higher co-dimension locus in the base. 

For the sheaf J 7 , one can find this singular locus by determining where the defining 
map g in (I4.18P goes to zero over the Calabi-Yau (I2.22p . In this case, there are 4 
point-like "instantons" , that is, points in the base for which J 7 jumps in rank. Such 
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The integrated error, t^ h , for SU (n) bundles on the Quartic K3, 
( 12.221) . Shown above are the results for 1) the stable SU{2) bundle 
V stable defined in (|3.17p . 2) the semi-stable SU(3) bundle ^ sem i- stable 
defined in (I4.9p . 3) the simple unstable sum of line bundles 0(— 1)© 
0{— 1) and 4) the unstable SU(2) bundle "V unstable defined in (|4.15p 
with a Harder- Narasimhan filtration exhibiting sheaf singularities. 
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instantons are described non-perturbatively in heterotic M-theory by a dissolved M5- 
brane on the K3 surface [63l 16^]. In homogeneous coordinates, these point^f] are 



(4.21) 



To explore the behavior of the connection around, say, po, we chose a geodesic P(s) 
that passes through this point. The approximate Calabi-Yau metric used for the 
geodesic was computed at k g = 15, and the parameter s is normalized such that it 
coincides with the path length. 

The geodesic is a solution of the differential equation 



Po 


= [i = 


-3.2279 


- 1.6903? 


-0.7510 


- 2.6031* 


: 1.3874 + 


2.8736*], 


Pi 


= [i = 


-3.2279 


+ 1.6903i 


-0.7510 


+ 2.6031* 


: 1.3874- 


2.8736*], 


P2 


= [i = 


-1.6951 


- 0.3844* 


1.6095 - 


0.5912* : 


-1.2183 + 


0.6535*], 


P3 


= [i = 


-1.6951 


+ 0.3844i 


1.6095 + 


0.5920*' : 


-1.2183 - 


0.6535*']. 



d 2p\. 



ds' 



lA dP»(s)dP v (s) 



ds 



ds 







(4.22) 



On a Kahler manifold, parallel transport does not mix holomorphic and anti-holomorphic 
coordinates. In other words, the Christoffel symbol 



V 1 



9 lh dj9kh 



(4.23) 



has no mixed-index components. Therefore, the equation for the geodesic simplifies 
to 



d 2 P l 



dPHs)dP k (s) 



ds 2 ' ~ JK ds 
We normalize the initial velocity vector as 



ds 







1 = 



dPHs)dPHs) 



ds ds 



s=0 



(4.24) 



(4.25) 



so that the geodesic parameter s equals the metric path length. To numerically solve 
for P(s), we used Newton's method with a step size of 10~ 3 . After each integration 
step, we projected back onto the K3 hypersurface. Some sample points along our 
chosen geodesic are 



P(0) = Po, 
P(0.05) = [1.005 + O.OOli 
P(l) = [l.640 + 0.372i 
P(2) = [2.127 + 0.968i 



-3.228 - 1.690* 
-3.228 - 1.690z 
-3.228 - 1.690z 



-0.749 - 2.593i : 1.387 + 2.879i] , 
-0.426 - 0.058i : 1.117 + 3.389i] , 
-0.528 + 2.7Ui : 0.608 + 3.50H] . 



(4.26) 



3 Such degeneracies can be easily found using computational algebraic geometry packages such as 

onus?]. 
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Figure 11: The eigenvalues of g lJ Fq along the geodesic P(s). The indecom- 
posable, unstable bundle "^unstable is the thick blue line, and one of 
its point-like instantons is located at s = 0. For comparison, the 
eigenvalues for the balanced connection on the sum of line bundles 
(I4.13D are given by the red line, while that for the stable bundle 
(I3.17P is shown by the green line. 



We now compute g^Fq for V unstab i c at each point of the geodesic. Recall that 
^unstable remains a bundle everywhere. We want to know: what happens to the color 
matrix g^Fq near the J-"-sheaf singular point P(0)? Since the vector bundles we 
are considering here are of rank 2, there are two real eigenvalues. After untwisting 
with the determinant line bundle, they are ±A(P(s)). We plot these eigenvalues in 



Figure llj Note, since V unstab i e is an SU(2) bundle, the two eigenvalues are A and 
—A. We plot only the positive eigenvalue. For comparison, we also compute the color 
matrix for the sum of line bundles O(-l) © 0(+l) and the stable bundle, V sta bi e 
from (13.17p along the same geodesic. Since we are interested in the accuracy of our 
numerical approximation to the slope-unstable connections, we list here the details of 
this preliminary calculation (in the sense that we have performed only a few iterations 
of the non- converging T-operator). 

• The Calabi-Yau metric (which determines the Kahler form u as wells as the 
geodesic P) was computed at degree k g = 15. 

• The metric T-operator was iterated 15 times, numerically integrating over 2,000,000 
points. 
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The Hermitian Yang-Mills connection was computed at degree kn = 7. 

The connection T-operator was iterated 30 times, numerically integrating over 
1,000,000 points. 



By inspecting |Figure lTj we see that the results for the stable bundle V stable and 
the unstable sum of line bundles 0(— 1) © 0(+l) behave smoothly along the geodesic. 
The eigenvalues for the stable bundle V stab i e (green dashed line) are close to zero, as 
they should be for a slope-stable bundle with ci(V stab i e ) = 0. Likewise, the eigenvalues 
for the reducible bundle 0(— 1) © 0(+l) are approximately constant and equal to the 
expected value 

i-A = Cl (0(l)) = l. (4.27) 

More interestingly, as predicted, the eigenvalues for the unstable bundle V unstab i e have 
a pole at -P(O), the location of the point-like instanton. Note that due to the log- 



scale in Figure 11, the width of this pole appears somewhat elongated. In reality, this 
divergence is very close to a 5-function, 5(s = 0). Away from that point, ^A ~ 1. 



This eigenvalue was precisely the value that was obtained in |Figure 10 by the general 



numerical integration (which did not include the points (I4.2ip ). 

As a final observation, note that the "dip" in the eigenvalues, where A — > before 
assuming the constant value, is a function of the fact that the bundle V unstable is 
not defined globally as a direct sum. That is, near the singularity at s = the 
eigenvalues of the color matrix have the form ( A _ A ) , whereas sufficiently far from the 
singularity, for s > 0.9, we find that the color matrix has "switched" the position of its 
positive and negative eigenvalues; that is, ( ~ 1 1 ) . In transitioning between these two 
configurations, the eigenvalues must go to zero. We conjecture that the failure of the 
eigenvalues to decompose globally into + and — values is in one-to-one correspondence 
with the fact that V unstab i e is not defined as a direct sum of line bundles but, rather, 
as an indecomposable rank 2 bundle. There is no physical significance to the point 
where these eigenvalues approach zero. In fact, we will see in the next subsections 
that as one increases the accuracy of our computation (specifically, higher iterations 
of the T-operator and of the ^//-twisting) the position of this zero approaches the 
5- function spike at s = 0. 



4.4.1 Iterating the T-Operator 

It is clear that the T-operator (13. 6 j) on the matrices H a p does not converge for V unstable . 
As noted in ISubsection 3.1\ the T-operator converges if and only if the bundle is 
Gieseker stable and V unstab i e is manifestly Gieseker unstable^. Indeed, since the T- 

4 To see that V unsta bic is not Gieseker stable, note that from (|3.7|) . Po(i) (J~)(n) — xO? 7 ® 0(n)) = 
2n(2 + n) while P 0{1) (V)(n) = x(V ® 0{n)) = 2(1 + n 2 ). Hence, for n > 1, P (i)(V) -< Po{i){F) 
and V is Gieseker unstable. 
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operator has to reproduce the field-strength singularity at the location of the point- 
like instanton, some components of the H a s matrix must grow beyond bound as we 
approximate the balanced connection. 

To see this behavior, consider higher iterations of the T-operator. We approxi- 
mate the field strength Fq at fixed degree = 5, and compute the T-operator by 
numerically integrating using 500,000 points. The results of the first 20 iterations 
of the T-operater are plotted in |Figure 12} Despite the fact that the T-operator is 



clearly failing to converge, and, hence, the results of ITheorem 41 do not apply, the 
numerical evidence presented here strongly suggests that the balanced connection is 
still a well-defined limit as we continue iterating the T-operator. 
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Figure 12: The eigenvalues A(P(s)) of g^ Fq along the geodesic P{s). The 
point-like instanton is located at s = 0. The Hermitian Yang-Mills 
connection was approximated at degree ku = 5, and the bundle 
T-operator was iterated between and 20 times. 



As a final exploration of this behavior, we perform the computation of the con- 
nection at the limit of our accuracy, that is, for both a high degree of twisting and a 
large number of iterations of the T-operator. 



4.4.2 Degree of Twisting 

In this subsection, we compute the color matrix g^F^j at degrees ranging from ku = 2 
to ku = 8. The connection T-operator was numerically integrated using 1,000,000 
points. For small degrees ku, some entries of H a p grow quickly at each iteration, and 
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the iteration must be stopped before the machine precision becomes insufficient. This 
is why we performed only 7 and 15 iterations at kn = 2 and kn — 3. For larger 
degrees kjj > 4, the matrix entries grow slowly enough that this is not an issue. The 



results are shown in Figure 13 
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Figure 13: The eigenvalues A(P(s)) of g K 1 Ffj along the geodesic P(s). The 
point-like instanton is located at s = 0. The different lines corre- 
spond to different degrees ku at which the Hermitian Yang-Mills 
connection was approximated. The divergent peak at s = is ap- 
proaching a 5-funtion 5(s = 0) as kn — >■ oo. 



For a bundle of the form ^m^tabie in (I4.15p . we conjecture that as ku — > oo the color 
matrix can be decomposed as sum of a divergent piece (due to the sheaf J 7 in (I4.17P 
in its Harder- Narasimhan filtration) plus a constant eigenvalue contribution (due to 
the "smoothing" of its Harder-Narasimhan filtration to O(-l) © 0{1) far from the 
points (jOTjl ). That is, 



(point-like instanton) + (constant A elsewhere). 



(4.28) 



This conjecture explains the convergence of t^ h to a constant value shown in Figure 10 
Specifically, that 

lim T kH ( Vunrtahle) = 1™ Tfc H fo(-l) © 0(+l)) = 2 (4.29) 



for regions away from the point-like instantons. 
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In summary, the results of the past few sections have given substantial evidence 
that even in the cases where a bundle V does not solve the Hermitian Yang-Mills 
equations (and, hence, the T-operator does not coverge), the generalized Donaldson 
algorithm still produces a physically relevant connection. In particular, in view of 
the results in Figure 13 it seems possible that even in the case of an unstable vector 
bundle for which the color matrix g % ^Fq diverges, the generalized Donaldson algorithm 
is providing an accurate approximation of the field strength Fx. In the limit as ku 
increases, as well as the iterations of the T-operator, we find that while the T-operator 
is not converging, the resulting connection not only correctly displays the field strength 
singularities due the sub-sheaf singularities, but also the smooth asymptotic behavior 
away from any singular loci. 

With these observations in hand, we conclude that it is possible to apply the 
generalized Donaldson algorithm to an arbitrary holomorphic vector bundle and to 
determine its slope- stability properties by inspection of the results, such as those 
shown in Figure 10 This observation provides a useful new tool for determining the 
existence of supersymmetric vacua in heterotic compactifications. 



4.5 Rank Three Bundles on the Quintic 

Having learned what behavior to expect from stable, semi-stable, and unstable bundles 
under the generalized Donaldson algorithm, we now present examples involving higher 
rank bundles on a Calabi-Yau threefold. This is the type of geometry that is ultimately 
of interest to us in M = 1 supersymmetric heterotic compactifications. 

Specifically, we consider the deformed Fermat Quintic, denoted P 4 [5], defined by 

u 5 + v 5 + x 5 + + z 5 + \uvxyz = , (4.30) 

where u, v, x, y, z are the homogeneous coordinates of P 4 . As was the case on the 
Quartic K3, the Kahler cone of the Quintic is 1-dimensional and all line bundles C 
on X can be written as 0(n) for some integer n. The global sections of 0(n) on X 
can be defined as a polynomial space whose dimension is given by 

{0 n< 

( n + 4 ) 0<n<5 . (4.31) 

(T)-(Y) n>5 

With these definitions in hand, one can compute a basis of polynomials of H°(X, V eg) 
C kH ) of the degree given by (I4.3ip . 

As in previous sections, we will investigate three types of vector bundle that are 
respectively 1) stable, 2) semi-stable and 3) unstable. We will consider two different 
types of unstable bundles- one an indecomposable rank 3 bundle and the other an 
unstable sum of three line bundles. The monad bundles are defined as follows. 
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4.5.1 A stable bundle 



The stable 577(3) bundle V sta bi e is given by 

/"stable 







0(-2 



i ©3 



■» -1 



i®6 



V, 



stable 



where the monad map is defined (via a map on sections) as 



stable 



5u+24v+3x+7y+18z 
9u+22v+l2x+14y+3z 
6x+3y 
10u+12j/+10^ 
13u+26v+18x+13y+12z 
\ 20u+28v+5x+12y+17z 



lln+14j;+18a;+14j/ 
8«+4z+9z 
3u+v+2x+13y+18z 

6v+llx+9y 
5u+llt;+6a:+9i/+12z 
12u+Uv+28x+6y+7z 



29z 12u+16v+llx+13y+l6z \ 
20u+33'u+9x+18j/+342 » 
24«+23^+46a;+17y+38^ 
36u+21d+122;+372;+28z 
3«+18x+15y+122 i 
10u+4?/+62 / 



(4.32) 



(4.33) 



As in the previous sections, since the Kahler cone of the Quintic is one-dimensional, 
it is straightforward to verify that V sta bi e is slope-stable by checking that H°(X, A V) = 
for k = 1, . . . 2 (Hoppe's criterion [54l [62]). Applying the generalized Donaldson al- 
gorithm to the computation of the connection on V sta bie, we find that the T-operator 
converges as expected and that the error measure, (13.251) . converges smoothly to zero. 



The results are shown in Figure 14 



4.5.2 A semi-stable bundle 

A semi-stable SU (3) bundle is defined by the following short exact sequence 

/semi-stable 



— > 0(-3) 
where / semi -stable IS 

/ 

/semi-stable 



■*0(-l) 







V 



semi-stable 







(4.34) 



Au 2 +uv+llux+15vx+8x 2 +lbvy+xy+l6uz+2lvz+26xz+14yz+36z 2 \ 
5u 2 +4uv+19ux+9vx+8uy+6vy+l0xy+y 2 +2Auz+21vz+17xz+9yz 
(2u 2 +6uv+20vv+l5ux+5vx+l4:uy+12vy+19xy+5vz+4z 2 ) 
(4« 3 + 10u 2 v+13uvx+15uvy+uxy+16vxy+uy 2 +5xy 2 +5y 3 
\ +6u 2 z+3uvz+6v 2 z+7uxz+llvxz+9uyz+8vyz+9xyz+hy 2 z+8vz 2 +Qxz 2 +Ayz 2 ) ) 



(4.35) 

In this case, the bundle admits a subsheaf, O C V scmi _ sta bi e with fi(0) = /i(V scmi _ s t a bic)- 
Hence, V semi _ sta bic is properly semi-stable, not polystable, and will not admit a solution 
to the Hermitian Yang-Mills equations. Its Jordan-Holder filtration ( 14. 5 j) is given 
simply by 

F®0, (4.36) 
where the rank 2 sheaf J 7 is defined by 



-> T -> C(l)® 3 -> 0(3) ^0 . 



(4.37) 



Because the S-equivalence class containing V scmi _ stab ie contains a poly-stable repre- 
sentative, this semi-stable bundle may be brought arbitrarily close to a solution of 
the Hermitian Yang-Mills equations by varying the bundle moduli (that is, bringing 
V S emi-stabie closer to the direct sum bundle, ( 14.36!) ). 
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4.5.3 Unstable bundles 



As in previous sections, we will compare the behavior of two unstable bundles. 

1. An unstable direct sum: 

First, consider the simple unstable direct sum defined by 

V sum = 0{-l) © 0(-l) © 0(2) . (4.38) 

V sum is a direct sum of stable objects, but because the three line bundles do not 
have the same slope, the sum is unstable. As in (14.7j) . the color matrix is given 
by 

9 fj F fj =y -1 j . (4.39) 

The error measure, (]3.25p . is predicted to lie at r^ H = 4. 

2. An indecomposable unstable bundle: 

We can define an unstable SU (3) bundle as follows, 

/unstable . 

— ► 0(-2) > 0(-l)® 4 © 0(1) — ► V unstablc — ► (4.40) 

with 

(x+6y+18z \ 
13ii+llx+9i/ \ 
7u+llv+4z _ (4-41) 

(4u 3 +10« 2 v+13m>:r+15irey+n£y+16iJxy+ui/ 2 +5:ry 2 +5y 3 +6u 2 ,z+3m;z / 
+6v 2 z+7uxz+llvxz+9uyz+8vyz+9xyz+5y 2 z+8vz 2 +6xz 2 +4yz 2 ) / 

Here V unsta bie is destabilized by the sub-bundle 0(1) with /i(C(l)) > 0. The 
graded sum, (I4.5| . associated with its Harder-Narasimhan filtration is given 
simply by 

£©0(1), (4.42) 
where Q is a rank 3 sheaf with c±(Q) = —1 defined by 

-)• 0(-2) -)• 0(-l)® 3 -> £ -> . (4.43) 

As discussed in Subsections 14. 31 and 14. 4[ we would not expect that since |ci(£?)| = 
|ci(0(l))| = 1, the error measure should approximate Tk H = 2. 

In these calculations, the metric on the Quintic threefold was computed at degree 
kg = 8, the metric T-operator was iterated 10 times and the numerical integration 
was carried out with 2,000,000 points. In addition, the connection T-operator was 
iterated 10 times, the numerical integration used 1,000,000 points and the Tk H integral 
was computed numerically with 100,000 points. 

The comparison between the stable, (I4.32p . semi-stable, (14.34 p . and the two un- 



stable bundles, ( I4.38P and f )4.40p . is given in|Figure 15 
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0123456789 

k H 

Figure 1J±: The integrated error measure t^ h on the Quintic threefold for the 
stable bundle defined in (|4.32p . 

5 Conclusions and future work 

Donaldson's algorithm [38 | |39 | 14*0] has been shown to be an important tool in the com- 
putation of Ricci-flat metrics on Calabi-Yau manifolds. Numerical implementations of 
this algorithm were given in [3SJ EHJ |50j H2] and a generalization of these techniques to 
Hermitian metrics on holomorphic vector bundles was proposed in (12} H3] . In this pa- 
per, we have presented a systematic and efficient method to implement the generalized 
Donaldson algorithm and to numerically compute Hermitian Yang-Mills connections 
satisfying (11. ip . We illustrated this by, first, computing the Ricci-flat Calabi-Yau 
metrics on the Quartic K3 surface and the Quintic threefold and, second, over these 
manifolds, calculating the Hermitian Yang-Mills connection on several holomorphic 
vector bundles defined via the monad construction [5T1 [521 ESI ES] • 

In addition to showing that the algorithm converges to the Hermitian Yang-Mills 
connection for a slope-stable bundle, we presented results demonstrating that even in 
the case of unstable vector bundles, for which the generalized T-operator (I3.6P does 
not converge, the algorithm produces a physically relevant connection. In Section HJ 
we showed that the connection produced by the generalized Donaldson algorithm in 
the case of an unstable bundle can be understood in terms of the Harder-Narasimhan 
filtration [51] of an unstable sheaf by semi-stable subsheaves. Furthermore, as shown 
in Subsection I4.4[ in cases where the nitrations contain sheaves with curvature singu- 
larities, the presence of such singularities can be accurately described by the algorithm. 
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Figure 15: Comparison of different bundles on the Quintic threefold. The sta- 
ble bundle, (|4.32j) . exhibits the behavior of a Hermitian Yang-Mills 
connection and the error t^ h decreases as the degree ku is in- 
creased. Meanwhile, the unstable, (|4.40p . and strictly semi-stable, 
(|4.34p . bundles approach a fixed, non-zero constant value. The re- 
sults for the semi-stable bundle depend on the choice of bundle mod- 
uli and can be made to converge to a constant arbitrarily close to 
zero. 
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These results may shed light on the mathematical study of non-Hermitian Yang-Mills 
connections [68] . 

Importantly, our results for unstable bundles allow us to apply the generalized 
Donaldson algorithm to arbitrary vector bundles arising in heterotic string compact- 
ifications and use it to determine whether such geometries admit supersymmetric 
vacua. The problem of deciding whether or not a given holomorphic vector bundle is 
slope-stable is a notoriously difficult one. Particularly challenging is the fact that the 
difficulty of a stability analysis generally increases rapidly with the dimension, h 1 ' 1 , 
of the Kahler cone. One of the great advantages of the algorithm presented in this 
paper is that, unlike in a standard analytic analysis, such numerical calculations can, 
in principle, be performed with essentially equal ease in arbitrary h ' x dimensions, 
thereby providing an important new tool in the analysis of supersymmetric heterotic 
vacua. 

In recent work, [561 [57], the dependence of the slope- stability property of a vector 
bundle on the Kahler and bundle moduli has been explored from the point of view 
of four-dimensional effective field theory. The presence of "stability walls" separating 
stable and unstable regions of Kahler moduli space has been shown to have a variety 
of interesting physical and mathematical consequences [SSJ EH]. In future work [70] . 
we will use the numerical techniques presented here to further explore the moduli 
dependence of supersymmetric solutions. 

Finally, it is worth noting that the results of this paper are an important step 
forward in the computation of observable quantities of particle physics, such as the 
matter-field Kahler potential and physical Yukawa couplings. We hope to explore 
these applications in the future. 
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